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Abstract 

Background: Heat and air pollution are both associated with increases in mortality. However, the interactive effect 
of temperature and air pollution on mortality remains unsettled. Similarly, the relationship between air pollution, air 
temperature, and social deprivation has never been explored. 

Methods: We used daily mortality data from 2004 to 2009, daily mean temperature variables and relative humidity, 
for Paris, France. Estimates of chronic exposure to air pollution and social deprivation at a small spatial scale were 
calculated and split into three strata. We developed a stratified Poisson regression models to assess daily 
temperature and mortality associations, and tested the heterogeneity of the regression coefficients of the 
different strata. Deaths due to ambient temperature were calculated from attributable fractions and mortality 
rates were estimated. 

Results: We found that chronic air pollution exposure and social deprivation are effect modifiers of the 
association between daily temperature and mortality. We found a potential interactive effect between 
social deprivation and chronic exposure with regards to air pollution in the mortality-temperature 
relationship. 

Conclusion: Our results may have implications in considering chronically polluted areas as vulnerable in heat 
action plans and in the long-term measures to reduce the burden of heat stress especially in the context of 
climate change. 
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Introduction 

Among environmental determinants of health, climatic 
effects are of high concern, especially given the growing 
body of literature concerning the impacts of climate 
change. Ambient temperature has long been recognized 
as a physical hazard, and is associated with a wide range 
of adverse health effects [1]. Increases in summer ambient 
temperatures are associated with increases in mortalities 
over a city specific threshold. Consequently, future heat- 
related mortality is likely to increase in the context of cli- 
mate change [2,3]. 
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Identification of factors that increase vulnerability to hot 
temperature in human populations has become a growing 
public health concern [4,5]. Some studies have explored 
different aspects of heat-related mortality or morbidity 
heterogeneity across different levels. These heterogeneity 
factors or effect modifiers in relation to temperature and 
mortality represent the heat vulnerability factors which 
can be observed at the individual or at the community 
level. At the scale of the individual, many studies have re- 
ported increased vulnerability of elderly individuals to hot 
temperatures [6-14]. The effects of gender on temperature 
impacts have seen mixed evidence [12,13]. At the commu- 
nity level, vulnerability of populations has been shown to 
be influenced by factors such as low density of green 
spaces [15], poor urban design and planning [16], urban 
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heat islands [17], and the modification effect by air pollut- 
ants [18,19]. Ozone in particular has been studied as an 
effect modifier in the mortality-temperature relationship 
[20,21], while the link between temperature and particu- 
late matter (PM) [22], sulfur oxides (SOx) [23] or Nitrogen 
oxides (NOx) [24] has been less explored. Furthermore, 
urban air pollution has never been studied from a chronic 
exposure perspective but only using daily levels. Some 
studies have also focused on community-based social 
characteristics that may influence mortality and morbidity 
outcomes through effects of community organization 
[25], social cohesion or social networks [26], and social 
deprivation using an aggregated index [27] . 

The interaction between ambient temperature and air 
pollution has been investigated in air pollution mortality 
time series studies and, to a lesser extent, in temperature 
mortality time series analyses [28]. To our knowledge, no 
data are available on the relationship between temperature 
and mortality effect modified by chronic air pollution ex- 
posure and socio-economic factors (or ecological social 
deprivation) assessed jointly. 

The objective of the study is to identify whether and 
how the magnitude of the effects of mean temperature 
on all-cause mortality were modified by chronic air pol- 
lution exposure (here nitrogen dioxide (N02) represent- 
ing urban traffic), social deprivation, and a combination 
of these two factors. For this purpose, we studied the 
city of Paris (France) where N02 long-term average con- 
centrations vary substantially across the city according 
to traffic density (i.e. the main source of N02 emissions), 
and where the different city neighbourhoods host popula- 
tions with contrasted socio-economic profiles. We hypoth- 
esized that a better understanding of these vulnerability 
factors should provide relevant information for developing 
public health programs targeting the most vulnerable pop- 
ulations and territories. 

Methods 

Study setting and small-area level 

Paris, the capital city of France, has a population of 
roughly 2.25 million inhabitants. The spatial scale used 
was the French census block (i.e. IRIS - a French acronym 
for "blocks for incorporating statistical information"), 
which constitute the smallest census unit areas in France, 
designed by the National Census Bureau (INSEE - Institut 
national de la statistique et des etudes economiques), for 
which aggregate data is available. The city of Paris is 
subdivided into 992 census blocks with a mean population 
of 2,199 inhabitants (range = 0-5,456 inhabitants) and a 
mean area of 0.11 km 2 (range = 0.009-5.4 km 2 ). 

Mortality and population data 

We considered all deaths occurred in the city of Paris (ex- 
cluding woods of Boulogne and Vincennes) for residents 



older than 35 years old from May to August for the years 
2004-2009 included. All-cause mortality data were pro- 
vided by the death registry of the city of Paris. Individual 
information on age, sex, date of death, and census block 
of residence was available for each case of death. For con- 
fidentiality issues it was not possible to distinguish causes 
of mortality, thus external causes of deaths could not be 
excluded. The analysis included only subjects older than 
35 years old at the time of death to minimize this bias 
because accidental causes of death are dominant in sub- 
jects under 35 years old [29]. We obtained the number of 
population in each stratum from the INSEE. 

Ethical approval was obtained from the French commis- 
sion on data privacy and public liberties (CNIL - Commis- 
sion Nationale Informatique et Liberte). 

Climate data 

Daily mean outdoor temperatures and relative humidity 
were obtained from Meteo-France, as measured at the 
Montsouris station in Paris and were computed using 
data of the corresponding period (summers from 2004 
to 2009). 

Air pollution data 

Annual N02 concentrations were modeled at a grid scale 
of 25x25m throughout the period 2004-2009 (only sum- 
mers) by the local association for the monitoring and 
the study of air quality (AirParif). We used a dispersion 
model (ESMERALDA) to produce annual N02 concen- 
trations at a fine spatial resolution. A description of the 
methods used to produce N02 levels at the French cen- 
sus block is provided in Additional file 1. Chronic N02 
exposure at the census block scale was defined as the 
average of annual N02 concentrations from 2004 to 2009. 
We did not include other air pollutants such as ozone due 
to the lack of data at this spatial scale. 

Community based socio-economic characteristics: social 
deprivation 

To characterize the socio-economic characteristics at a 
community level, we computed an aggregated, multidi- 
mensional deprivation index, fitted at the census block 
scale [30]. We provided a summary description of the 
deprivation index and its categorization in Additional 
file 1. 

Statistical analyses 

Social deprivation was defined and stratified according 
to the 3-classes deprivation index described above. N02 
chronic exposure was categorized into three groups ac- 
cording to the terciles of its distribution. To create a 
double stratification with sufficient statistical power, we 
stratified the chronic N02 exposure into two strata ac- 
cording to the median of N02 chronic concentrations (i.e. 
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we obtained a total of 6 strata for this double stratifica- 
tion). We also assessed age-related vulnerability (i.e. num- 
ber of daily deaths for the > 65 years people compared to 
the number of daily deaths for the <65 years) and sex- 
related vulnerability (number of daily deaths for men ver- 
sus the number of daily deaths women). 

First, in a crude model, we developed stratified Poisson 
regression models to assess daily temperature-related 
mortality associations with daily mean temperatures and 
daily death counts [31]. We used cubic B-splines (5 
knots) of time to control for secular trends in the mor- 
tality series [9]. Seasonal patterns of mortality were con- 
trolled by including a quadratic function represented by 
the day of the season (1 to 123) [32]. Sensitivity analyses 
were performed by using a cubic polynomial for season- 
ality. We used natural cubic splines (2 knots) to con- 
sider the non-linearity of the temperature-mortality 
relationship. Daily levels of humidity were incorporated 
into regression models as possible confounding vari- 
ables based on evidence in the literature [1,33,34]. Then, 
we developed stratified models for each category [35]. 
The threshold for statistical significance was set at p < 
0.05 and all the tests were two-sided. Model assump- 
tions and validity were verified graphically (i.e. with 
quintile-quintile and partial autocorrelation function 
plots). A white noise test [32] was also used to ensure 
that no auto-correlation remained in the residuals. Esti- 
mates of exposure (daily temperature) - response (mor- 
tality) functions - were expressed as Relative Risks (RR) 
and their confidence intervals (CI) at 95% for each 
temperature degree. Relative Risks for the relation be- 
tween mortality and temperature were estimated relative 
to the daily mean number of deaths for the entire period 
[34]. We assessed the heterogeneity of temperature-mor- 
tality associations across different strata using the method 
developed by Payton et al. [36] to compare the regressions' 
coefficients between different strata. This test was con- 
ducted on all regression coefficients across the strata. 



from attributable fractions (AF), using heat-related mor- 
tality relationships (RR) described in the previous section 
[37,38]. We considered only days with daily mean 
temperature related to a RR strictly greater than 1 for any 
temperature value (corresponding to 280 days in total). 
AF were calculated for each temperature degree (Ti°: 
represents a temperature unit) and using the equation 
[(RR(Ti°)-l)/RR(Ti 0 )]. 

The Attributable Number of deaths (AN) for the period 
2004-2009 (months of May, June, July, and August), for a 
given temperature variable, was then estimated from the 
following equation: 

max 

AN = ^2 A¥ ( T i) x MDC x ND(r,-) 

i—min 

In this equation, i is the temperature degree Celsius, 
AF(r,) is attributable fraction for temperature degree 
(Ti), MDC is the mean observed daily death count, and 
ND(T,) is the number of days for which temperature was 
i (°C). We then divided the total number of attributable 
deaths for the period 2004-2009 by 6 to obtain an aver- 
age attributable number of deaths by summer. We used 
the RR estimated within the strata of interest to estimate 
the AR in those strata. We calculated confidence intervals 
at 95% of attributable number of deaths. Then we com- 
pared the number of deaths attributable to temperature in 
each of the groups stratified by chronic air pollution ex- 
posure (3strata), social deprivation (3 strata) and the two 
simultaneously (6 strata). We made sure that there was no 
homogeneity in the temperature-mortality associations 
across different strata, as described in the last section. Fi- 
nally we calculated attributable mortality rates, by dividing 
the total number of deaths attributable to temperature in 
each stratum (by summer) by the total population in each 
of these strata. The final results were mapped to visualize 
heat vulnerability according to social deprivation and 
chronic air pollution exposure. 



Calculation of deaths attributable to temperature 

We calculated deaths attributable to ambient tempera- 
tures for each stratum. Attributable deaths were calculated 



Results 

46,056 deaths were registered during the study period, 
and the mean age at death was 68 years (SD = 10.21). 



Table 1 Summary statistics for health outcomes, climate and air pollutants data (Paris, 2004-2009) 


Variable 


Mean 


Minimum 


25th Percentile 


Median 


75th Percentile 


Maximum 


Standard deviation 


Daily death count (n) 


33.3 


3 


27 


33 


42 


93 


4.5 


Minimum Temperature (°Q* 


12.5 


0 


9 


13 


16 


21 


4.2 


Mean Temperature (°Q* 


17.3 


/ 


14 


17 


20 


26 


4.4 


Maximum Temperature (°Q* 


20.9 


8 


1/ 


21 


24 


35 


5.2 


Relative Humidity (%)* 


68.4 


37 


61 


69 


76 


94 


10.7 


Daily N02 (ug/m 3 ) 


48.9 


16 


38 


47 


58 


121 


14.3 


Chronic N02 (ug/m 3 ) 


52.5 


38.7 


47.88 


51.8 


56.1 


81.1 


6.8 



*Summary statistics for weather variables (minimum, mean and maximum temperatures as well as relative humidity) are based on daily data. 
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As presented in Table 1, daily death count ranged from 
3 to 93 (mean, 33.3, IQR: 27 to 42). Mean temperatures 
ranged from 7°C to 26°C (mean: 17.3°C, IQR: 14°C to 
20°C), while atmospheric pollution assessed by daily 
N02 ranged from 16 to 121 ug/m3 (mean: 48.9 ug/m 3 , 
IQR: 38 ug/m 3 to 58 ug/m 3 ) (Table 1). Chronic N02 
exposure ranged from 39 ug/m 3 to 81 ug/m 3 (mean: 
52.5 ug/m 3 , IQR: 48 ug/m 3 to 56 ug/m 3 ). From the 992 
census blocks in our study, we had 54 census blocks 
with missing data concerning social deprivation which 
corresponded to 160 deaths. These census blocks are non- 
residential (activity and miscellaneous) and with few resi- 
dents. Descriptive statistics of the social deprivation index 
and chronic air pollution are presented in (Additional 
file 1: Table SI). As shown in Additional file 1: Figure SI, 
chronic N02 concentrations were slighdy higher in the 
most favoured census blocks. In the crude model (for the 
whole population), mortality was significantly associated 
with daily mean temperature (RR for each temperature 
unit and for all strata are shown in Additional file 1: Table 
S2 and Table S3). The association between temperature 
and mortality was U-shaped. Applying the central point 
estimate of the temperature-mortality relationships to the 
observed mean temperatures, an average number of 121 
attributable deaths by summer (May to August) was esti- 
mated. Deaths attributable to mean temperature are pre- 
sented with their 95% confidence intervals in Additional 
file 1: Table S4. 

As presented in the Figure 1, we observed a gradient 
in the raw deaths attributable to temperature according 
to social deprivation: the higher the degree of social 
deprivation, the higher the number of deaths attributable 
to mean temperature. Likewise, higher chronic air pollu- 
tion exposure as assessed by N02 levels was also identi- 
fied as a vulnerability factor in the relationship between 
temperature and mortality. We observed a significant 
difference in regression coefficients (between strata) 
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Figure 1 Number of deaths attributable to mean temperature 
by social deprivation (3 strata) and chronic air pollution 
exposure (3 strata). 



separately according to social deprivation and chronic 
air pollution (p < 0.05). 

Finally, we calculated mortality rates (Table 2), including 
the population living in each strata to consider the differ- 
ences in population density across strata. We found some 
heterogeneity (p < 0.10) between ambient temperature 



Table 2 Mortality rates attributable to summer 
temperatures (per 1 00 000) by strata 

Strata Population Mortality rates p for 

{per 100 000 f heterogeneity* 



Total 
Age 

Under 65 years 

More than 65 years 

Sex 

Female 

Male 

Social Deprivation 

Low Social Deprivation 

Medium Social 
Deprivation 

High Social Deprivation 

Chronic air pollution 
exposure 

Low Chronic N02 
Exposure 

< 50.6 ug/m 3 

Medium Chronic N02 
Exposure 

50.6-55.8 ug/m 3 

High Chronic N02 
Exposure 

> 55.8 ug/m 3 

Double stratification 
Low Chronic Exposure 
Group 

Low Social Deprivation 

Medium Social 
Deprivation 

High Social Deprivation 

Double stratification 
High Chronic Exposure 
Group 

Low Social Deprivation 

Medium Social 
Deprivation 

High Social Deprivation 



2234105 5.37 [5.01 ;5.73] 

1921330 0.78 [0.62;0.88] 

312774 3357 [31.6535. 81] 

1161734 5.51 [4.99;5.85] 

1072370 5.22 [4.76;5.67] 

692572 4.33 [3.31 ;5.23] 

781936 4.60 [3.595-41] 

743956 7.26 [6.74;7.85] 

737254 4.75 [4.1 3;5.22] 

795341 5.97 [5.36;6.32] 

670231 7.89 [7.28;8.14] 



370567 3.78 [2.87;5.03] 
381190 4.19 [3.345.25] 



375341 



6.92 [5.1 1 ;8.1 2] 



389872 3.59 [2.295.07] 
391432 5.36 [4.22;6.41] 

325701 9.82 [7.79;10.93] 



0.001 



0.51 



0.08 



0.03 



0.14 



0.07 



a For social deprivation and chronic air pollution strata, mortality rates were 
standardized for age and sex. 

b Homogeneity test of Chi- square Pearson, for central estimates. 
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and mortality according to age, social deprivation and 
chronic air pollution exposure. Sex was not an effect 
modifier in the relationship between heat and temperature 
(p = 0.51). For the double stratification, we found that, in 
the low chronic exposure group, social deprivation did not 
significantly modify the relation between heat and mor- 
tality (p = 0.14), while we found that social deprivation 
may have modified this relation in the low chronic ex- 
posure group (p = 0.07). We also conducted heterogen- 
eity tests across the two chronic exposure groups. We 
found a heterogenic effect (data not shown) only be- 
tween the strata of low chronic exposure/high social 
deprivation and the strata of high chronic exposure/ 
high social deprivation. 

Figure 2 presents the spatial distribution of social 
deprivation and chronic air pollution exposure. The 
spatial variations of both social characteristics are re- 
lated to deaths attributable to temperature as presented 
in Table 2. 



Discussion 

In this study, we assessed the effects of mean temperature 
and chronic air pollution (with N02 representing urban 
traffic) on all-cause mortality with additional attention 
paid to the interactive effects of social deprivation. To do 
this we conducted stratified time-series analyses. We 
found that chronic air pollution exposure modifies the 
association between daily temperature and mortality. 
We also found that social deprivation is a heat-related 
vulnerability factor. Finally we found that there is a poten- 
tial combined modification effect of social deprivation and 
chronic exposure to N02 with regards to heat-related 
mortality. 

Our results for the temperature-mortality relationship 
are comparable to other studies conducted in Paris [8,32]. 
Concerning social vulnerability to heat-related mortality, 
our results are consistent with other studies, using time 
series analyses [5] or using different methods, such as a 
case crossover design [11]. 
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Figure 2 Spatial distribution of social deprivation and chronic air pollution. 
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Ecological social vulnerability can be explained by an 
accumulation of different types of vulnerabilities for in- 
dividuals within a census block (e.g., social isolation [39], 
material conditions [40], poor urban design including 
micro heat islands [17] or poorer health [14]). In the same 
way, social mobility may also partly explain the observed 
social gradient, similar to previous studies looking at 
socio-economic scales [41]. Thus, populations with low 
socio-economic levels will be inclined to stay in their 
neighbourhoods, with a low daily mobility, accumulat- 
ing other heat vulnerability factors and with other con- 
sequences on certain determinants of health such as 
physical activity. 

We found that air pollution modifies the relation be- 
tween ambient temperature and daily mortality as found 
in previous studies using daily exposure [21-23,28,42]. It 
is possible that chronic N02 levels modify the effects of 
temperature on mortality. A range of studies have shown 
that N02 is consistently associated with many health 
outcomes such as cardiovascular effects [43,44], which 
can explain how the N02 chronic effects can make 
some populations more vulnerable to heat impacts. 

Our results suggest that ambient temperature has a 
greater impact on daily mortality on chronically polluted 
areas, especially if these areas are socially deprived. How- 
ever our results only suggest this double effect modifier 
and further studies may confirm this result. 

Methods for estimating the heat related mortality 
relationships inevitably rely on some assumptions. In 
our study, we chose to not take into account lag effects 
to simplify the interpretation of deaths attributable to 
temperature. We considered that contribution of lag ef- 
fects compared to the day of death was not differential 
according to the different strata characteristics. We did 
not take into account harvesting effect or mortality dis- 
placement. "We also did not consider intra urban varia- 
tions of temperatures. Future studies should consider 
these aspects. We considered the social deprivation at 
the ecological level, and we did not explored social in- 
equalities in mortality at the individual level, because of 
the lack of data. The use of a synthetic index rather than 
using independent ecological socio-economic measures 
allows us to represent an accumulation of social and ma- 
terial disadvantages [45] and these kinds of indexes are 
more useful to lead population interventions [27]. The 
division of neighbourhoods into census blocks aims to 
maximize their homogeneity in terms of population size, 
socioeconomic characteristics, land use and zoning. In 
Paris, the area of a block census is quite small so the 
ecologic bias is likely to be small. Finally, we did not 
consider the spatial auto-correlation in our analysis which 
is a result of the correlation between the spatial unit and 
the adjacent geographical areas. Variables observed at 
small area level are interdependent because proximity and 



linkages between neighbouring areas. This issue could 
influence the effect of social deprivation and chronic air 
pollution [46], and further studies could consider it by 
conducting spatial analyses such as a hierarchical Bayesian 
modeling approach. 

Conclusion 

In this study, we showed that areas which are chronically 
polluted by air pollution can be characterized as vulner- 
able to heat in the sense that they have more deaths attrib- 
utable to heat than less exposed areas. We also presented 
a potential combined vulnerability of social deprivation 
and chronic exposure to air pollution in the mortality- 
temperature relationship. Our results may have important 
implications considering chronically polluted areas as vul- 
nerable in heat actions plans (especially including adapted 
surveillance and warning systems) and in the long-term 
measures to reduce the burden of heat stress (as building 
regulations, urban planning or land-use changes), espe- 
cially in the context of climate change. However, further 
studies are necessary to determine whether similar results 
could be found in other settings. 
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